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q , We study the anomalous cross section lineshape of e e — ¥ DD with an effective field theory. 

■ Near threshold, most of the DD pairs are from the decay of ^>(3770). Taking into account that 

the non-resonance background is dominated by the ip(2S) transition, the produced DD pair can 
undergo final state interactions before it is detected. We propose an effective field theory for the 
low energy DD interactions to describe this final state interactions, and find that the anomalous 
lineshape of the DD cross section observed by the BESII collaboration can be well described. 
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" As the first charmonium state above the DD threshold, the resonance "0(3770) is different from other 

charmonia with lower masses. Since the 0(3770) decay into the open charm DD is allowed by the 
OZI rule, this dominant decay mode leads to a broad width up to 27.2 ± 1.0 MeV Obviously, the 
direct production process of e + e~ — > 0(3770) — > DD is useful for the study of the property of ■0(3770). 
In Ref. [2], BESII collaboration reported an anomalous behaviour of the cross section lineshape at the 
0(3770) mass region in e + e~ — > DD that cannot be described by a simple Breit- Winger of 0(3770). Such 
an observation has inspired interesting theoretical discussions |3j-|6[. In particular, it was found that the 
1 interfering effect between 0(3770) and ^^S 1 ) plays a very important role in understanding the anomalous 

lineshape of DD at the 0(3770) resonance [3, |4|. Such an interference can be recognized by a relative 
phase factor e 1 ^ which is introduced between these two resonances, and the phase angle turns out to 
be large in order to describe the anomalous DD lineshape. 

In principle, the phase factor e 1 ^ can come from the final state interactions of DD. Thus, it should 
O s \ be interesting to study the DD anomalous lineshape using the effective field theory to describe the DD 

final state interactions. This forms our motivation for this work. Near threshold, the DD pair produced 
in e + e~ — > DD comes from the decay of the 0(3770) and other non-resonance background processes. 
Once the DD pair was produced, it could undergo final state interactions before it converts into the final 
observed DD state. This could explain the relative phase between the "0(3770) and other non-"0(377O) 
amplitude, and provide a description of the DD lineshape. We note that there are several cases in which 
the final state interactions play important roles in the understanding of the cross section lineshapes [7l-fl0| . 

It is well-known that effective field theory is a useful tool to study the low energy hadron interactions. 
An effective field theory makes use of the Tailor expansion on the small ratio between the typical small 
scale p and the cutoff scale A. For example, in Chiral Perturbation Theory (ChPT), p is the momentum 
of the low energy pion or pion mass while A = M p ( 770 ) sets the cutoff scale of this effective theory. 
An effective field theory for the low energy DD is different from that for the low energy 7T7t interaction 
because the 0(3770) should be included explicitly into the effective Lagrangian. Besides the three-vector 
momentum of the D(L>) meson, another small scale 6 = M^( 3 77 ) — 2Mr> ~ 40 MeV also appears in the 
effective theory. This additional small scale will make the power counting to be different from that in 
ChPT. Systematic development of the effective field theory with resonances as intermediate states is still 
under exploration, and interesting discussions on this subject can be found in Refs. [iH fl2T]. 

In this work, we use the effective field theory to study the DD interaction in order to understand the 
dynamic details of the anomalous cross section lineshape observed by the BESII Collaboration [2|. 

At the beginning, we assume that the production of DD in e + e~ annihilation can be approximated by 
the vector meson dominance (VMD). It means that the cross section for e + e~ — > DD is dominated by 
intermediate vector meson productions via e + e~ — > 7* — > IZi — > DD, where IZi denotes any vector meson 
with isospin / = or I = 1. However, it is impossible to sum the contributions from all of Ri in reality. 
As a reasonable approximation, one can include the contributions from the vector mesons in the vicinity 
of the considered energy region, but neglect those far off-shell vector mesons. In the energy region of 



the BES data from 3.74 GeV to 3.8 GeV, one can expect that 0(3770) plays the most important role 
among all of Ri, while the contributions from all the other Ri can be treated as background. As shown 
in Ref. [3j, the contribution from ip(2S) dominates the background, while contributions from other states 
are negligible. Therefore, we only include the contributions from the resonances 0(3770) and %p(2S), and 
neglect those from the other resonances. Namely, ip(2S) would be the main background near the DD 
threshold. 

InVMD [HEil, 

the coupling between the vector meson and virtual photon can be described as 

pM 2 

C Vl = ^V^\ (1) 
Jv 

where is the vector meson field, is the photon field, and My is the mass of the vector meson. 
Setting the electron mass m e « 0, the coupling can be obtained as 
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where T ee is the electron-position decay width of V^, and a = 1/137 is the fine-structure constant. 

Once the DD pair is produced from the decay of the vector meson 0(3770) or 0(25), it can undergo final 
state interactions through the rescattering processes DD — > DD —>•••—» DD which can be described by 
the effective field theory. In the energy region which is concerned, the three-vector momentum p of the 
D(D) meson is small. Thus, it is possible to construct an effective field theory for the low energy DD 
interactions making use of the expansion on the small momentum p. Since the mass of 0(3770) is just 
above the threshold of DD, we need to include -0(3770) explicitly in the formulation. Near threshold, 
the D(D) meson can be treated as nonrelativistic. Thus, the interaction Lagrangian for the DD system 
with the quantum number J PC = 1 can be constructed as 



CjDD = ^g^ D D{D^D-VDW}-xP + ig^ DD {D^VD-VD^D}-xP, 
\ddY = fx{D^VD-VD^D}-{VD^D-D^VD] 

+/ 3 {D t r ! V5 - V.DVD} • {VL>VL> - DVV£>} + • • • , 
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with D = l D+ j, D=l"_ y (3) 

Here ip is the field operator of -0(3770), D(D') annihilates a D(D) meson while D'(D) crates a D(D) 
meson, t % is the Pauli matrix, and the ellipsis denotes other contact terms with more derivatives which are 
higher order terms. The first term in C^j^y accounts for the interaction in the isospin singlet channel, 
and the second term for the isospin triplet channel. The contributions from other resonances which are 
not included in the Lagrangian can be saturated into the contact term £(£>£>)2. Therefore, we take the 
coefficients such as /i , /3 to be complex, where the imaginary parts of them come from the width of the 
saturated resonances and the DD annihilation effect. With isospin symmetry, we only have to consider 
the terms for isospin singlet channel in C^ D ^y to study the DD final state interactions, since the DD 
pair comes from the decay of 0>(377O) and 0(25) in our approach. 

Now we come to the discussion on the power counting of this effective field theory. The tree level 
diagrams for the DD elastic scattering are shown in Fig. [T] Near the DD threshold, the denominator of 
the 0(3770) propagator can be expressed as 
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FIG. 1: Tree diagrams for DD ->■ DD. (a) DD ->■ V>(3770) Dfl, (b) contact interaction. 



where p is the magnitude of the three- vector momentum of the D(D) meson in the overall center of mass 
frame, T^ on - DD denotes the non-DD decay width of ?/>(3770), and M$ is the mass of ?/>(3770). The DD 
decay width of "0(3770) will be included through the summation of the D meson loops in the following. 
Because V>(3770) is close to the threshold of DD, we expect that P(i/j) is a t 0(jT 2 ). Taking the momentum 
power of the ipDD vertex into account, we find that Fig.[IJa) is at O(jo ). From the naive power counting, 
the leading contact terms have two derivatives, hence they are at 0(jp 2 ). However, in this naive power 
counting, we have assumed that the coefficients of the contact terms, i.e. fx, fy, ■ ■ ■ , are at order of G(p°). 
In some cases, especially when there are bound states or resonances near threshold, the coefficients of the 
contact terms can be enhanced. For example, in NN interaction the S-wave contact terms Cs scale as 
0(p -1 )[l5|. Another example is the NN interaction in 3 Po, where the leading contact term Czp can be 
promoted to C(p _2 )[16j. It is interesting to study whether the same enhancement mechanism takes place 
in the DD interactions since the resonance ^(3770) locates near the DD threshold. If fx is promoted to 
0(p~ 2 ) as Cs Po in NN interaction, then the corresponding tree diagram shown in Fig. [IJb) is at O(p ), 
which is the same as Fig. QJa). However, since we do not know the power of fx at the beginning, we 
then assume that the leading contributions to DD elastic scattering come from both Fig. [Ua) and (b). 
We will use the experiment data to determine fx and see whether it is promoted to 0(p~ 2 ). In this way, 
the DD scattering amplitude in the specific channel (J PC = 1 ,/ = 0) can be obtained by summing 
the ladder diagrams as shown in Fig. [2] which is equivalent to solve the Lippmann-Schwinger equation 
T = V + f VGT with the DD potential truncated at leading order. 




FIG. 2: The ladder diagrams for the DD interactions where the potential is truncated at the leading order. 
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Figure [3] illustrates the final state interactions between the produced DD, and loop integrals that we 
will encounter in Fig. [3] can generally be reduced to 
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where E = p 2 /Mo is the total kinematic energy of the DD system. We can then write down the amplitude 
for e + e~ — > DD as 



iM = iM a + iMb- 
To be more specific, the amplitude for process Fig. [SJa) reads 
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where 7 is the Dirac gamma matrix, k\ and hi are the incoming momentum of electron and positron, 
respectively, while p\ and p 2 are the outgoing momentum of D and D, respectively. G^p cannot be 
simply interpreted as the width of ^(3770), since it is a complex number. If /1 = 0, then G,p = = 
T non-DD + {Wd0 \3 /Md0 + \p D+ \yM D+ )gl DB /(l2nM^. 
The amplitude for Fig. EJb) can be written as 
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where M^^s) is the mass of ip{2S) which can be read from PDG [17|, and f^(2S) can be extracted by 
Eq. © using r v , (2S) ^ e+e - = 2.35 keV 

To proceed, we denote the cross section for e + e~ — > DD as cr B (s) which does not include the initial 
state radiation (ISR) effect. In reality, for a given energy -^/s, the actual cm. energy for the e + e~ 
annihilation is 



: yj s(l — x) due to the ISR effect, where xE^ eam is the total energy of the emitted 

photons. To order a 2 radiative correction in the e + e~ annihilation, the observed cross section cr ^ s at 
BESII can be related to our result a B through [18[ 



a ohs (s) = (1 + S VP ) 
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(11) 



where (1 + Svp) = 1.047, and the function f(x, s) is given by 
f(x,s) = (3x^ 1 S v+s + S H , 
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FIG. 3: The Feynman diagrams for e + e — > DD in our approach. 



Before fitting the BESII data with Eq. (TTT]) . we first discuss our treatment of r non_D ' D . It seems 
impossible to determine Y^ on ~ DD definitely in our fitting since Y 1 ^ 011 ' 00 is always accompanied by fx 
in our formula, and any change of r non_I313 can be compensated by the tuning of fx . The experimental 
results on non-DD branching ratio of -0(3770) decay are still controversial fl9T - l2lj ]. In contrast, the 
next-to-leading-order (NLO) pQCD calculation expects the non-DD decay branching ratio to be at most 
about 5% 22], while an effective Lagrangian approach estimates that the D meson loop rescatterings 



into non-DD light vector and pseudoscalar mesons leads to about 1% non-DD branching ratios [23[. A 



similar calculation by Ref. 24j also confirms such a nonperturbative phenomenon. Another experimental 
observation is that so far most of the well-measured non-DD decay modes of -0(3770) are found rather 
small, namely, their branching ratios are either at the order of 10~ 3 — 10~ 4 , or only an upper limit is 
set fV% , 

Taking all these facts into account and for the purpose of studying the dominant DD channel, we set 
pnon-DD £ Q zerQ j n Qur fitting as a leading approximation. We have checked that the fitting results 

are almost unchanged even though we set the non-DD breaching ratio of -0(3770) decay to be at the 
order of several percents. Our fitting results for the DD cross section lineshape are presented in Fig. [4] 
The fitted parameters are 

M$ = 3.7674(44) GeV, = 4.25(13) x 10~ 3 , 

9ipDD = 15.43 ±2.7, 94,(2S)dd = -6.9 ± 3.6, 

fx = 2059.2 ± 534 (GeV -2 ). (13) 

The fitting quality is x 2 /d.o.f. = 23.44/22 ps 1.07. The imaginary part of fx from the fitting is almost 
zero, which is consistent with the fact that the DD annihilation will be highly suppressed in the heavy 
quark limit. It seems that the magnitude of fx is large, and it is worth noting that this value is at the 
order of (Md/S) 2 , where 5 = Af^( 37 7o) — 2Mb- Therefore, the large fx is consistent with our previous 
assumption that fx will be promoted by the low energy scale <5, up to (D(p~ 2 ). With the fitting parameters, 



5 



we can obtain 



1 



(M 2 - AMlo) 3/2 + (M 2 - AM 2 D+ f/' 2 



27.58 ±11.07 MeV, 



r PP = ^aM^(-) 2 = 165.6 ± 10.1 eV. 

3 J-0 



(14) 



The values of M^ and T ee we obtain are somewhat smaller than those given by PDG2010 |l7j, while 
is consistent with the PDG value [l7j- We note that the 0(3770) resonanceparameters, such as M^, 
and r ee , are still under improvement as shown by the PDG2012 revisions 

In Fig. El we also present the Born cross section a B (s) for e + e~ — > DD which is denoted by the solid 
curve. The dashed and dotted lines are for the neutral and charged D meson pair Born cross sections, 
respectively. Combining the results shown in Figs. 0] and [5l we find that the anomalous cross section 
lineshape could be originated from the interferences from the 0(25) pole and DD final state interactions. 
Since the 0(25) pole is relatively isolated due to its relatively narrow width in comparison with the 
mass gap between 0(25) and 0(3770), the relative phase between the 0(25) and 0(3770) amplitudes is 
likely to be produced by the DD final state interactions. Although our calculation cannot determine the 
absolute value for the possible non-DD decay branching ratio of ^(3770), it is constructive to recognize 
the important role played by the final state DD interactions that cause the deviation of the e + e~ — s- DD 
cross section at the -0(3770) mass region from a Breit-Wigner shape. This analysis is useful for our 
further understanding of the "0(3770) non-DD decays as a manifestation of possible nonperturbative 
QCD mechanisms. 
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FIG. 4: The observed cross sections for e + e~ — > D°D° (left plot) and e + e~ — ¥ D + D~ (right plot). The solid line 
is the fitting result in our approach shown in FigO the dashed line shows the contribution from Fig(3](a), and 
the dotted line shows the contribution from FigEJ (b) . The data are from BES0]. 



In summary, we have proposed an effective field theory for low energy DD interactions, in which we 
have included the resonance "0(3770) and an additional small scale 6. It is found that the small scale 5 
will promote the coefficient of the contact term j\ to be 0(p~ 2 ). So the leading DD interaction potential 
in this specific channel would come from the S'-channel 0(3770) exchange and the contact term f\. With 
the leading DD potential, we then sum the ladder diagrams to describe the DD final state interaction 
as shown in Fig. [31 We find that we can describe the anomalous cross section lineshape of e + e~ — > DD 
observed by the BESII Collaboration Q using the effective field theory. This approach should be useful 
for our further understanding of the 0>(377O) non-DD decays which could be sharing the same dynamic 
origin as the DD cross section lineshape anomaly as emphasized in Refs. (23l. [25|. 
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FIG. 5: The Born cross section cr s (s) for e + e — > DD. The solid line is for e + e — > DD, the dashed line is for 
e+e" -s> D°D°, and the dotted line is for e+e - 
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